Introduction to the Bethe Ansatz I 
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A few years after the formulation of quantum mechan- 
ics, Heisenberg and Diractl discovered that one of its won- 
drous consequences was the key to the age-old mystery of 
ferromagnetism. They found that the laws of quantum 
mechanics imply the existence of an effective interaction, 
Jij Si ■ Sj , between electron spins on neighboring atoms 
with overlapping orbital wave functions. The exchange 
interaction, as it has become known, is caused by the 
combined effect of the Coulomb repulsion and the Pauli 
exclusion principle. This spin interaction was soon recog- 
nized to be the key to a microscopic theory of ferromag- 
netism and many other cooperative phenomena involving 
electron spins. _ 

In 1931 Hans BetheB presented a method for obtain- 
ing the exact eigenvalues and eigenvectors of the one- 
dimensional (ID) spin-1/2 Heisenberg model, a linear 
array of electrons with uniform exchange interaction be- 
tween nearest neighbors. Bethe's parametrization of the 
eigenvectors, the Bethe ansatz, has become influential to 
an extent not imagined at the time. Today, many other 
quantum many body systems are known to be solvable 
by some variant of the Bethe ansatz, and the method 
has been generalized and expanded far beyond the ad 
hoc calculational tool it was originally. Unlike the sim- 
ulation of a classical model system, most computational 
approaches to quantum many-body systems require a fair 
amount of analytical work up front. This requirement is 
true for Monte Carlo calculations, renormalization group 
approaches, the recursion method, large-scale numerical 
diagonalizations, and for the Bethe ansatz when used in 
a computational context. 

In spite of its proven importance and wide range of 
applications, the Bethe ansatz is rarely discussed in text- 
books on quantum mechanics and statistical mechanics 
except at the advanced level. The goal of this column is 
to introduce the Bethe ansatz at an elementary level. In 
future columns we plan to discuss some of the extensions, 
generalizations, and applications of the method, to which 
many workers have made important contributions. Our 
emphasis will be on computational applications for which 
the method is less well known than for its usefulness as 
an analytical tool. 

The Bethe ansatz is an exact method for the calcu- 
lation of eigenvalues and eigenvectors of a limited but 
select class of quantum many-body model systems. Al- 
though the eigenvalues and eigenvectors for a finite sys- 
tem may be obtained with less effort from a brute force 
numerical diagonalization, the Bethe ansatz offers two 
important advantages: (i) all eigenstates are character- 



ized by a set of quantum numbers which can be used to 
distinguish them according to specific physical proper- 
ties; (ii) in many cases the eigenvalues and the physical 
properties derived from them can be evaluated in the 
thermodynamic limit. 

The Hamiltonian of the Heisenberg model of spins 
S n = (S%,S%,S%) with quantum number s = 1/2 on a 
ID lattice of N sites with periodic boundary conditions 
Sjv+i = Si is given by 

N 

H = — J S„ • S„ + i 
1 



n=l 
N 



(1) 



where S„ = S% ± iS® are spin flip operators. H acts on 
a Hilbert space of dimension 2 N spanned by the orthog- 
onal basis vectors \ai . . . ctn), where a n =t represents 
an up spin and a n —I a down spin at site n. The spin 
commutation relations (with H = 1) are 

[Sn>S n i] — ±S^5 nn i, [S^,S n /] — 2S^5 nn i . (2) 

The application of the operators , on a vector 
|ai . . . ojv) yields the results summarized in Table B. We 
can use these rules to express H as a real and symmetric 
2 N x 2 N matrix whose eigenvectors can be computed by 
standard diagonalization algorithmsa. From the eigen- 
vectors, the physical quantities of interest can be calcu- 
lated by evaluating the expectation values (matrix ele- 
ments) for the corresponding operators. 

TABLE I. Rules governing the application of the spin op- 
erators on the basis vectors \cri . . . on) with a n =|, J.. 
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The Hamiltonian matrix can be written in block di- 
agonal form if we perform, prior to the numerical di- 
agonalization, one or several basis transformations from 
{|(Ti . . . ctat)} to a symmetry- adapted basis. These trans- 
formations reduce the computational effort needed for 
the remaining numerical diagonalization and make it pos- 
sible to handle larger system sizes. 

The Bethe ansatz is a basis transformation that does 
not have to be supplemented by a numerical diagonaliza- 
tion. In principle and usually in practice, this alternative 
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procedure removes the cap on system sizes. However, the 
implementation of the Bethe ansatz entails calculational 
challenges of its own. Depending on the specifics of the 
application, they can be met by analytical or computa- 
tional methods as we shall see in the following. 

For the Heisenberg model, two symmetries are essential 
for the application of the Bethe ansatz. The rotational 
symmetry about the z-axis in spin space, which we have 
chosen to be the quantization axis, implies that the z- 
component of the total spin Sf. = Yl n =i IS conserved: 
[H, Sj,] = 0. According to the rules of Table |, the oper- 
ation of H on | £7j . . . un) yields a linear combination of 
basis vectors, each of which has the same number of down 
spins. Hence, sorting the basis vectors according to the 
quantum number Sj- = N/2 — r, where r is the number 
of down spins, is all that is required to block diagonalizc 
the Hamiltonian matrix. 

The block with r — (all spins up) consists of a single 
vector \F) = | f ... f). It is an eigenstate, H\F) — 
E \F), with energy E = -JN/A. The N basis vectors 
in the invariant subspace with r = 1 (one down spin) are 
labeled by the position of the flipped spin: 



= S-\F) n=l,. 



(3) 



To diagonalize the r = 1 block of H, which has size 
N x N, we take into account the translational symmetry, 
i.e., the invariance of H with respect to discrete transla- 
tions by any number of lattice spacings. Translationally 
invariant basis vectors can be constructed from the vec- 
tors in ([}]) by writing 



N 

E 



e ikn \n) 



(4) 



for wave numbers k = 2Trm/N, m = 0, . . . , N — 1. (The 
lattice spacing has been set equal to unity.) The vec- 
tors are eigenvectors of the translation operator with 
eigenvalues e lk and are also eigenvectors of H with eigen- 
values 



E - E = J(l - cosfc), 



(5) 



as can be verified by inspection (see Problem 1). The vec- 
tors (^) represent magnon excitations, in which the com- 
plete spin alignment of the ferromagnetic ground state 
\F) is periodically disturbed by a spin wave with wave- 
length A = 2n/k (see Problem 2). 

In the invariant subspaces with 2 < r < N/2, the 
translationally invariant basis does not completely diag- 
onalizc the Hamiltonian matrix even if we take into ac- 
count further symmetries of (|l|) such as the full rotational 
symmetry in spin space or the reflection symmetry on the 
lattice. Here the Bethe ansatz is a powerful alternative. 

For r — 1, the case we have already solved, we now 
proceed somewhat differently. Any eigenvector in the 



r = 1 subspace is a superposition of the basis vectors 
(I): 



N 
n=i 



(0) 



The eigenvector is a solution of the eigenvalue equa- 
tion H\ip) — E\ip) if the coefficients a(n) satisfy the linear 
equations 

2[E - E }a{n) = J[2a(n) - a{n - 1) - a(n + 1)] (7) 

for n — 1,2,... , N and with periodic boundary condi- 
tions a(n + N) = a(n). N linearly independent solutions 
of (0) are 



a{n) 



ikn 



k = — m, 

N 1 



0,1,... ,N-1. 



(8) 



If we substitute the coefficients a(n) into (g), we obtain 
(after normalization) the magnon states (0) with energies 

§■ 

The distinctive features of the Bethe ansatz begin to 
emerge when we apply the same procedure to the case 
r = 2. The task is to determine the coefficients a(n\, n%) 
for all eigenstates 

\ip) = ^2 a(ni,n 2 )\ni,n 2 ), (9) 

l<ni<n 2 <JV 

where \n\,n2) = S~ S~ \F) are the basis vectors in this 
subspace of dimension N(N — l)/2. Bethe's preliminary 
ansatz for the coefficients is 

a( ni ,n 2 ) = Ae^ kini+k2n ^ + A 1 e t{kin2+k2ni) . (10) 

We might wish to set A — A', use the same values for 
k\,k 2 as in (||), and interpret the wave function as a 
superposition of two magnons. However, the result would 
be an overcomplete set of N(N + 1)/2 nonorthogonal and 
nonstationary states. Superimposed spin waves are in 
conflict with the requirement that the two flipped spins 
must be at different sites. The eigenvalue equation for 
(^) translates into N(N — l)/2 equations for as many 
coefficients a(ni,n 2 ): 

2[E - E ]a(m,n 2 ) = J\Aa(n\, n 2 ) — a{n\ — 1, n 2 ) 
— a(rai + l, n 2 ) - a(m,n 2 — l) - a(m,n 2 + l)] 
for n 2 > ni + 1, (H a ) 

2[E — E ]a(ni,n 2 ) — J[2a(ni,n 2 ) — a(ni— 1, n 2 ) 
— a(m,n 2 +l)] 

for n 2 =ni + l. (Hb) 

Equations (11a) are satisfied by a(ni,n 2 ) in ( |l(i| ) with 
arbitrary A, A' , k±, k 2 for n 2 > fix + 1 and for n 2 = fi\ + 1 
provided the energy depends on fci, k 2 as follows: 



E-En 



J (1 — coskj). 

.7=1,2 



(12) 
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Equations ( lib ), which are not automatically satisfied by 
the ansatz (|10|), are then equivalent to the N conditions 

2a(ni, ri\ + 1) = a(n\, rii) + a(n\ + 1, ni + 1) (13) 



obtained by subtracting ( lib ) from ( 11a ) for n 2 = ni + 1. 
In other words, ajni, n 2 ) are solutions of Eqs. jll| ) if they 
have the form ( |l0| ) and satisfy (|ll|). The conditions jl3| ) 
require a modification of the amplitude ratio, 



.4 



e i(k 1 +k 2 ) -i_ 1 _ 2e ikl 

e i(k 1 +k 2 ) _|_ 1 _ 2e i?C2 



(14) 



This requirement is incorporated into the Bethe ansatz 
as extra phase factors 

0(711,713) ^ e i(kini+k2n2+ ^ 012) + e l{kin2+k2ni + ^ 821 \ (15) 

where the phase angle £q 2 = —O21 = depends on the as 
yet undetermined k%,k2 via ( |l4| ) or, in real form, via 



9 ki 
2 cot - = cot — 
2 2 



cot ■ 



(16) 



The quantities k\ , &2 will henceforth be referred to as 
momenta of the Bethe ansatz wave function (^|) with co- 
efficients (|l5|). 

Two additional relations between k\,k2, and 9 fol- 
low from the requirement that the wave function be 
translationally invariant, which implies that 0(711,712) = 
a(n,2,ni + N). This condition is satisfied by the coeffi- 
cients @ if e tklN = e ie , e lk2N = e~ 10 . Equivalently, we 
can write (after taking logarithms) 



Nkx = 27rAi + f 



Nki = 2ttA 



9, (17) 
1} are called Bethe 



where the integers A 2 ; G {0, 1, . 
quantum numbers. 

The remaining task is to find all (Ai,A2) pairs which 
yield solutions of Eqs. ( |l6| ) and (p7[), known as the Bethe 
ansatz equations. Every eigenstatc in the r = 2 subspace 
can be found in this way. For any solution k\, A^, 9, the 
(non- normalized) eigenvector has coefficients of the form 
( |l5| ) . The expression (JlJ) for the energy and the relation 



k = k\ 



A 2 ) 



(18) 



for the wave number k are reminiscent of two superim- 
posed magnons. The magnon interaction is reflected in 
the phase shift 9 and in the deviation of the momenta 
k\ , &2 from the values of the one-magnon wave numbers 
as given in (||) . We shall see that the magnons either scat- 
ter off each other or form bound states. Note that the 
momenta k\ , specify the Bethe ansatz wave function 
(^) , while the wave number k is the quantum number as- 
sociated with the translational symmetry of H and exists 
independently of the Bethe ansatz. 

The analysis of the complete r = 2 spectrum will 
demonstrate the usefulness of the Bethe quantum num- 
bers for distinguishing eigenstates with different physical 



properties. The allowed (Ai,A2) pairs are restricted to 
0<Ai<A2<A r — 1. Switching Ai and A2 simply 
interchanges k\ and and produces the same solution. 
There are N(N + l)/2 pairs that meet this restriction, 
but only N(N— 1)/2 of them yield a solution of Eqs. jl^ ) 
and (0). The solutions can be determined analytically 
or computationally. Some of them have real k\,k2 1 and 
others yield complex conjugate momenta, &2 =k\. 

We first find all the solutions and then interpret them. 
The (Ai,A2) pairs which yield solutions for N = 32 are 
shown in Fig. 1. We begin with the class C\ of states 



Xi + A 2 = 3N/-2 



A, + A, = N/2 



























Ai + A 2 = N 
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FIG. 1. The allowed pairs of Bethe quantum numbers 
(Ai, A2) that characterize the ^(A^ — l)/2 eigenstates in the 
r — 2 subspace for N = 32. The states of class C\,Ci, and 
C3 are colored red, white, and blue, respectively. 

for which one of the Bethe quantum numbers is zero, 
Ai = 0, A 2 = 0, 1, . . . , N — 1. There exists a real solution 
for all N combinations, fci = 0, &2 = 2tt\2/N, 9 = (see 
Problem 3). 

Next consider the class C2 of states with nonzero Ai , A2 
which differ by two or more: A 2 — Ai > 2. There are 
N(N — 5)/2 + 3 such pairs. All of them yield a solu- 
tion with real k\,k2- To determine these solutions, we 
combine Eqs. (fl6|), (17), and ( fi8| ) into a single nonlinear 
equation for k\\ 

Nk\ k\ k — ki , . 

2 cot - = cot— -cot -. (19) 

2 2 2 y ' 

The solution k\ of (^|) for a given wave number k = 
(2-Kn/N), n = 0, 1, ... ,7V - 1 is substituted into @ 
and yields &2 and 9. Some of the solutions can be found 
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analytically (Problem 4) , and others must be determined 
numerically (Problem 5). 

The remaining class C3 of states has nonzero Bethe 
quantum numbers Ai, A 2 which either are equal or differ 
by unity. There exist 27V — 3 such pairs, but we will 
see that only N — 3 pairs yield solutions of (|l6| ) and (17). 
Most of the class C3 solutions are complex. To find them, 
we write 



k 

fci = - + IV, 



k 2 = - - iv, 



IX, 



(20) 



and use Eqs. ( |i6| ) and ([17]) for fixed k to obtain the rela- 
tion 



c, 

Ci 

Cj . x,=x, 

C, , X,=X,-l 



cos — sinh(A r w) = sinh[(7V — l)v] + cos <f> sinh v, (21) 

where <fi = 7r(Ai — A2), and \ — Nv is inferred from the 
solution. It is sufficient to consider v > 0. The energy 
( |l2| ) of any complex solution is rewritten in the form 



E - E n = 2J j 1 — cos — coshw 



(22) 



Inspection of ( pl| ) shows (see Problem 6) that a complex 
solution exists for A2 = Ai if 



Ax + A 2 = 2,4,... , --2,^+2, 



,2JV-2. 



(23) 



For A2 = Ai + 1 a complex solution exists if 



A!+A 2 = A,A + 2,... ,iV/2-l,3/V/2 + 3, ... ,2jV-A+2, 

(24) 



where A w \/N /it. In the latter case, there exist addi- 
tional real solutions of ([n]) if 



Ai + A 2 = 3, 5, . . . , A - 2, 2N - A + 2, . 



,2N 



(25) 



The (Ai, A2) pairs in (P3|)-(P5|) account for N— 4 solutions 
in class C3. There is one more class C3 solution. This 
state, which has k — tt and Ai = A 2 = N/4, is easily 
missed in the numerical analysis, because fci, k% have an 
infinite imaginary part (see Problem 7) . 

The complete r = 2 excitation spectrum (E — Eq)/J 
versus k for a system with N = 32 spins as obtained by 
the analytical and computational procedures outlined in 
Problems 3-7 is plotted in Fig. 2. The N states of class 
C\ form a branch with exactly the same dispersion rela- 
tion (||) as the magnon states that populate the r = 1 
subspace. This degeneracy is a consequence of the con- 
servation of the total spin St (see Problem 8). The class 
C 2 states are spread in a regular pattern over a region in 
(k, i?)-space. They are nearly free superpositions of two 
one-magnon states. 

The excitations belonging to classes C% and C 2 can 
be characterized as two-magnon scattering states. The 



FIG. 2. Excitation energy (E — Eq)/J versus wave number 
k of all N(N — l)/2 eigenstates in the invariant subspace with 
r — 2 overturned spins for a system with N = 32. States of 
class C\ are denoted by red circles, states of class C2 by open 
black circles, and states of class C3 by blue squares if A2 = Ai, 
or blue diamonds if A2 = Ai + 1. 



effect of the magnon interaction on these states is visu- 
alized in Fig. 3. It shows all N(N - l)/2 - N+ 3 class 
Ci and class C 2 states for N = 32 in comparison with 
the N(N + l)/2 two-magnon superpositions, where the 
momenta kj, j = 1, 2 in (|12| ) are replaced by one-magnon 
wave numbers kj = 2Ttrrij/N, rrij = 0,1,... , JV — 1. The 
magnon interaction manifests itself as a modified excita- 
tion energy of the two-magnon scattering states. Note 
that the interaction energy approaches zero when either 
ki or fc 2 goes to zero. The class Cj states can then be 
interpreted as exact superpositions of a ki = magnon 
and a /« 2 7^ magnon. 

As N increases, the energy correction due to the 
magnon interaction diminishes for the class C 2 states as 
well and vanishes in the limit N — ► 00. The two-magnon 
scattering states then form a continuum with boundaries 
E — Eq = 2J(1 ± cosfe/2), which coincide with those of 
the continuum of free two-magnon states. 

For the class C3 states, the effects of the magnon inter- 
action are much more prominent, and they do not disap- 
pear in the limit N — > 00. In the (k, -E)-plane of Fig. 2, 
these states lie on a single branch with dispersion (see 
Problem 9) 



E — E = —(1 — cosfc) 



(26) 



below the continuum of two-magnon scattering states. 
They are the two-magnon bound states. 

The bound state character of the class C3 states man- 
ifests itself in the enhanced probability that the two 
flipped spins are on neighboring sites of the lattice. 
This property of the wave function is best captured in 
the weight distribution |a(rti,rt 2 )| of basis vectors with 
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FIG. 3. Excitation energy (E — Eo)/J versus wave number 
k of all two-magnon scattering states (classes Ci and C2 from 
Fig. 2) for a system with N = 32 in comparison with the 
noninteracting magnon pairs (+). 

flipped spins at sites n\ and ni. In Fig. 4 we plot 
\a{ni, n2)| versus 712 — n\ for a sequence of class C3 states 
between A; = and k — n. The distribution is peaked at 
n<i — ni = 1. Its width is controlle d by the imaginary 
parts of k\,k2,9 in the coefficients (fl5|). The smallest 
width is observed in the bound state at fc = tt, that is, 
for Ai = A2 = N/4, whose fci,fc2 have an infinite imag- 
inary part (see Problem 7). In this case, all coefficients 
in (^|) with 7J2 7^ ni + 1 are zero, which implies that the 
two down spins are tightly bound together and have the 
largest binding energy. For the adjacent bound state with 
quantum numbers Ai = A2 = N /4 — 1, the exponential 
dependence of the weight \a(ni,n,2)\ on the separation 
ri2 — n% can be worked out analytically (see Problem 9). 

The width of the distribution |a(ni, 77,2) | increases as fc 
decreases, and the binding of the two down spins loosens. 
For finite N, the exponential factors disappear when the 
distribution has acquired a certain width, and the Bethe 
ansatz solutions switch from complex to real. In con- 
trast, the distribution |o(ni,?i2)| for scattering states is 
always broad and oscillates wildly in general. However, 
for some combinations of fci, fc2, a smooth distribution en- 
sues, which has its maximum when the two down spins 
are farthest apart (ri2 = n\ + N/2) (see Problem 10). 

The formation of bound states and scattering states of 
specific elementary excitations exist in a large variety of 
physical contexts. But only in rare cases such as this one 
can the nature and the properties of such compound ex- 
citations be investigated (analytically and computation- 
ally) on the level of detail made possible by the Bethe 
ansatz. 

We now present the Bethe ansatz for an unrestricted 
number r of overturned spins. We generalize (0) and 



FIG. 4. Weight distribution |a(m,7i2)| versus distance 
ri'z — n\ of the two down spins of class C3 states at 
k = (2ir/N)n, n = 4, 8, . . . , N/2 for TV" = 128. 

expand the eigenstates in the form 

\tp) = ^2 a(ni,... ,n r )\ni,... ,n r ). (27) 

l<ni<...<n,.<JV 

The subspace has dimension Nl/[(N — r)!r!]. The gener- 
alization of ([l5|) for the coefficients in terms of r momenta 
kj, and one phase angle 6ij — —Qji for each (fci, k-i) pair 
is as follows: 

a(m, . . . , n r ) = exp I i ^ k Vj nj + - 8 ViVj . 

ves r \ 3=1 Kj ) 

(28) 

The sum V G S r is over all r! permutations of the labels 
{1,2, ... ,r}. For r = 2 the two permutations are the 
identity (1, 2) and the transposition (2, 1), which produce 
the two terms of (|l5|). The consistency equations for 
the coefficients a(ni, . . . , n r ) inferred from the eigenvalue 
equation H\ijj) = E\ip) are 

r 

2[E-E ]a(m,... ,n r ) = X! [ a ( n i>>-- > n r) 

i=l a=±l 

-a(ni, . . . ,rii + a, . . . ,n r )}, (29a) 
for rij+i > rij + 1, j = 1, . . . , r, and 
2[E - E ]a(m, ... ,n r ) = 

r 

J ^ ^ [a(ni,... ,n r )-a(ni, . . . ,m+a,... ,n r )\ 
+jy^[2a(ni,... ,n r )-o(Tii,... ,71^-1,71^1,... ,n r ) 

a 

-a(ni, . . . ,nj a ,rij a+ i + 1, . . . ,7i r )] (29b) 
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for rij a+ i = rij a + 1, rij+i > n 3 - + 1, j ^ j a . The coeffi- 
cients a(ni, . . . , n r ) are solutions of (29) for the energy 



TABLE II. Bethe ansatz solutions for N = 4, r = 2. 



£-£ = J^(l-cosfcj) (30) 
i=i 

if they have the form (^ ) and satisfy the conditions 
2a(n l5 . . ;Tij a ,Tij a +l, ■ ■ .,n r ) = 

a(ni, . . ■,nj a ,nj a , . . .,n r )+a(ni, . . ., n ]a +1, n jQ +1, . . .,rt r ) 

(31) 

for a = 1, . . . , r just as we have shown for r = 2. These 
conditions relate every phase angle 9ij to the (as yet un- 
determined) kj in (p7|): 



e i(fc«+fe 3 ) + 1 _ 2e ik * 

g to ij = 

e i{k t +kj) _|_ i _ 2 e *fej 

An equivalent relation in real form reads 

y 2 



(32) 



(33) 



The translational invariance of ( |27| ) implies that the 
coefficients fl2"s| ) satisfy the relation a(ni,... , n r ) = 
a(n2, ... , n r ,ni + AT). Consequently, we must have 



2 cot = cot — cot — , i, j = 1, . 



r I i 

3=1 «3 
r 

+ E kv(j-i) n j + k-p/ r (ni + N), (34) 

3=2 

where the relation between the permutations on the left 
and the right is V'(i -l)=Vi, i = 2, . . . , r; V'r = VI. 
If we take into account that all terms not involving the 
index V'r — VI cancel, we are left with r additional 
relations between the phase angles and the momenta: 



3#* 



%3 j 



i=l, 



(35) 



where Ai S {0, 1, . . . , A^ — 1} as in (|i7j). What remains to 
be done is to find those sets of Bethe quantum numbers 
(Ai,... ,A r ) which yield (real or complex) solutions of 
the Bethe ansatz equations ( [3 3 ]) and ( |35| ) . Every solution 
represents an eigenvector ([27|) with energy ( j30|) and wave 
number 



2tt 
~N 



(36) 



i=l 



The complete set of Bethe ansatz solutions for systems 
with A^ = 4, 5, and 6 spins can be read off Tables 0, 
HI, and |v|, respectively. The solutions are for the in- 
variant subspace with r — N/2 (5f. = 0) for even A^ or 
r = (N- 1)/2 (Sf. = 1 /2) for odd N. For given momenta 



S T 


AiA 2 


2k/ ty 


fci 


k 2 


E - E 


2 

















1 


1 


1 





tt/2 


1 


1 


2 


2 





7T 


2 


1 


3 


3 





3tt/2 


1 





1 3 





2tt/3 


4tt/3 


3 





1 1 


2 


7r/2 + ioo 


tv/2 — ioo 


1 



TABLE III. Bethe ansatz solutions for N = 5, r = 2. 



S T 


Ai A2 


5k/2iv 


fci 


k 2 


-E — Eo 


5/2 

















3/2 


1 


1 





2tt/5 


0.690983 


3/2 


2 


2 





4tt/5 


1.809016 


3/2 


3 


3 





6tt/5 


1.809016 


3/2 


4 


4 





8tt/5 


0.690983 


1/2 


1 3 


4 


1.705325 


3.321222 


3.118033 


1/2 


1 4 





tt/2 


3tt/2 


2 


1/2 


2 4 


1 


2.961962 


4.577859 


3.118033 


1/2 


1 1 


2 


2tt/5 + £1. 198913 


2tt/5 - il. 198913 


0.881966 


1/2 


4 4 


3 


8ir/5 + il. 198913 


87r/5 - il. 198913 


0.881966 



ki, the phase angles 9ij can be determined from (|32j) and 
the eigenvectors (p7|) by the substitution of fej , into 
(^8|). These solutions yield all the energy levels. Some of 
the levels are degenerate, namely those with one or sev- 
eral Bethe quantum numbers equal to zero. The remain- 
ing eigenvectors of any such level belong to the same St 
multiplet but have different r, i.e., different S^- All Bethe 
ansatz solutions for r < A^/2 (even N) or r < (N — l)/2 
(odd AT) can be inferred from the ones listed in Tables 0, 
III, IV by removing the momenta ki — one at a time. 
Each ki — removed from a solution in the r subspace 
yields a solution in the r + 1 subspace with the remaining 
ki unchanged. 

For large AT, the classification of the Bethe ansatz solu- 
tions becomes more and more intricate as r increases to- 
ward N/2, and finding them all becomes increasingly te- 
dious. A more modest, but nevertheless highly promising 
and useful goal in many applications, is to find selected 
solutions for very large systems (N — > 00) - solutions that 
determine specific physical properties (ground state en- 
ergies, magnetization curves, susceptibilities, excitation 
spectra) of the underlying model system. This approach, 
to which many authors have made important contribu- 
tions since 1931, will be explored in future columns. 

We conclude with a brief discussion of some Bethe 
ansatz solutions for r > 2 that are of particular im- 
portance in the context of the Heisenberg ferromagnet 
(J > 0). Earlier we have found that two magnons may 
form a bound state with a considerable binding energy 
(see Fig. 2). We have seen that in these states the 
two down spins are much more likely to be on nearest- 
neighbor sites than is the case for two-magnon scattering 
states. The fact is that three or more magnons can form 
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TABLE IV. Bethe ansatz solutions for N = 6, r = 3. 



S T 


Ai A2A3 


3k/n 


fei 


k 2 


k 3 


E — Eq 


3 




















2 


1 


1 








71-/3 


1/2 


2 


2 


2 








2tt/3 


3/2 


2 


3 


3 








7T 


2 


2 


4 


4 








4tt/3 


3/2 


2 


5 


5 








5tt/3 


1/2 


1 


1 3 


4 





1.419506 


2.769283 


2.780775 


1 


14 


5 





1.340040 


3.895947 


5/2 


1 


1 5 








27r/5 


8tt/5 


1.381966 


1 


2 4 








471-/5 


67r/5 


3.618033 


1 


2 5 


1 





2.387237 


4.943144 


5/2 


1 


3 5 


2 





3.513901 


4.863679 


2.780775 


1 


1 1 


2 





7t/3 + 10.732857 


tt/3 - i0.732857 


0.719223 


1 


5 5 


4 





5tt/3 + i0.732857 


5tt/3 - i0.732857 


0.719223 


1 


12 


3 





7r/2 + ioo 


7r/2 — ioo 


1 





3 


3 


il.087070 


-il.087070 


7T 


0.697224 





114 





tv/2 + ioo 


7r/2 — ioo 


7T 


3 





1 1 5 


1 


1.338006 + il.471688 


1.338006 - il.471688 


4.654369 


2 





1 5 5 


5 


1.628815 


4.945179 + il.471688 


4.945179 - il.471688 


2 





1 3 5 


3 


1.722768 


7T 


4.560416 


4.302775 



bound spin complexes with even lower excitation ener- 
gies. Specifically, for the subspace with r down spins, 
it can be shown that the lowest excited state at fixed 
wave number k is represented by the wave function with 
(complex) momenta given by 



k 1 k 
cot -f = r cot — 
2 2 



i(r-2 J + l) + 0(e- s > N ), (37) 



with Sj > 0, j — 1,... ,r. This solution generalizes 
the case r = 2 discussed in Problem 9 and yields exact 
solutions for N — ► 00. The dispersion of the resulting 
bound state branch with r < N/2 is 



E - E = ~(1 - cosfc). 



(38) 



The Bethe quantum numbers of any such state are char- 
acterized by \Xi — Aj+i| = 0, 1 for every pair of complex 
conjugate momenta ki+i = k* in (|3^). 

In part II of this series, the focus will be on the ID 
s = 1/2 Heisenberg antiferromagnet. We shall employ 
the Bethe ansatz to determine the ground state of this 
model, the spectrum of low lying excitations, and the 
calculation of transition rates for dynamical quantities. 



Suggested problems for further study 

1. The translation operator T shifts the local spin con- 
figuration to the left by one site with a wrap around 
at the ends, e.g., T|TTU) = |TUT>- Show that the 
states (Q) are eigenvectors of T with eigenvalues e lk 
and eigenvectors of H with eigenvalues (||). 

2. Each spin wave state contains one flipped spin. In a 
traveling wave it is located at each lattice site with 
equal probability. Show that the periodic nature 
of the disturbance in the spin alignment of a spin 
wave state is reflected in a 1 /TV-correction of the 
spin correlation function: 



IS, -S 



The nearest-neighbor correlation function defines 
an effective angle 6 between nearest-neighbor spins: 
(^|Sz • Si + x\ip) = icos(9. From @ with n = 1 
we see that the smaller the wavelength, the larger 
8 and the larger the energy (|5|) of the state. For 
k = the spins remain fully aligned, and the state 
is degenerate with \F). 

3. Show that the Bethe ansatz equations for Ai = 
0, A2 = 0, 1, . . . , N — 1 are solved by k\ = 0, &2 = 
2n\ 2 /N,e = 0. Use @ instead of which is 
singular in this case. 

4. Express (|l|) for k = in the form cot(iVfci/2) = 
cot(fci/2), which yields k x = -k 2 = 2lir/(N - 1) 
with integer I. For k = it, express ( [l9| ) in the form 
cot(A r fci/2) = cot(fci), which yields fci = tt — fc 2 = 
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21-k I (N—2). For any such solution k-y, k 2 found, use 
( |l2|) to determine the excitation energy and ( |l6| ) 
and ( |i~7| ) to determine Ax,A2- 

Solving (|l^) is equivalent to finding the zeros of the 
function f(x) = 2cot(-/Vx) — cot a; + cot(/c/2 — x), 
where x = k\/2. Standard subroutines such as can 
be found in Numerical Recipesu ask for an interval 
[a;i,X2] that contains exactly one zero of f(x) such 
that f{x\) ^ 0, f{x2) §s 0. To fine tune this pro- 
cedure, it is useful to preview the distribution of 
zeros by plotting f{x) for < x < tt and various 
N. 



6. Solving (|2lJ) is equivalent to finding the zero at v > 
(if one exists) of the function f(v) = k sinh(7Vw) — 
sinh([7V— l]v) — 6 sinhw, where k — cos(fc/2), S = 1 
for A2 = Ax, and S = — 1 for A2 = Ax + 1. Consider 
the case A2 = Ax. Show that a zero of f(v) exists if 
Nk > and N(k - 1) < 0, which implies ([||. A 
similar argument yields the allowed Bethe quantum 
numbers for the case A2 = Ax + 1. Find the zero 
of f(v) numerically for all allowed combinations of 
A2 = Ax and A2 = Ax + 1 and determine the energies 
via M). 



7. Insert (|20|) with k = tt, cp = 0, \ = vN into QlSj), 
and show that a(nx,rt2) = i(— l) ni ^n 2 -m +1 a s v ~ * 
00. Show that these coefficients satisfy ( |llb| ) with 
E — Eq = J. Show that the same coefficients result 
from the solution for Ax = A2 = 37V/4. This is the 
only r = 2 state whose Bethe quantum nximbers 
are not unique. 



8. Show that the operator 



=1 S„) corn- 
Show that 



mutes with H and Sf. by using 
S%\F) = S T {S T + l)\F) with S T = N/2. Conser- 
vation of the total spin St implies that all eigen- 
states are at least (2 St + l)-fold degenerate. The 
level with energy Eq has St — N/2, and hence its 
degeneracy is N + 1 . It is represented in all r sub- 
spaces. For r = it is the state \F), for r = 1 it 
is the one-magnon state with k = 0, and for r = 2 
it is the class C\ state with k — 0. The remain- 
ing one-magnon states in the r = 1 subspace have 
S T = N/2 - 1, which implies an (N - l)-fold de- 
generacy. They are represented in all r subspaces 
with 1 < r < N — 1. For r = 2, they are the class 
C\ states with k ^ 0. All other r = 2 states have 
S T = N/2 - 2. 

9. Show that for N — * 00, the solution of f(v) = 
with f(v) as defined in Problem 6 converges 
toward v(k) = — lncos(fc/2), which substituted 
into yields @. Show that cotfc^ 2 /2 = 

2cotfc/2 =F i. For k = tt — 4n/N, this solution 
is extremely accurate even for small N. Show 
that \a{ni,n 2 )\ oc cosh[(a; — l/2)AlnA], where 
< x — (n 2 - ni)/N < 1. Plot this distribu- 
tion versus x for several N and with appropriate 



normalization. Demonstrate that all coefficients 
a(nx, 712) with n 2 — n\ > 1 go to zero as N — > 00. 

10. Show that Eqs. @ and @ for A x ,2 = A/2=Fl are 
solved by fcx, 2 = 7r[l=Fl/(JV-l)] with (E-E )/J = 
4 sin 2 tt /{N — 1). This two-magnon scattering state 
is the highest r = 2 excitation energy. The wave 
function of this state has coefficients a(nx,^2) = 
2(_i)«i+n 2 s i n [ 7r (i/2 + m - n 2 )/(N - 1)]. Plot 
|a(nx, n2)| versus 712 — 77x and compare this weight 
distribution with that of a two-magnon bound 
state. 
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